Cumulants on the Wiener Space 



by Ivan NourdirEl and Giovanni Peccatlll 

Universite Paris VI and Universite Paris Quest 



Abstract: We combine infinite- dimensional integration by parts procedures with a recursive rela- 
tion on moments (reminiscent of a formula by Barbour (1986)), and deduce explicit expressions for 
cumulants of functionals of a general Gaussian field. These findings yield a compact formula for 
0^ ■ cumulants on a fixed Wiener chaos, virtually replacing the usual "graph/diagram computations" 
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■ 1 Introduction 

' The integration by parts formula of Malliavin calculus, combining derivative operators 

and anticipative integrals into a flexible tool for computing and assessing mathematical 
expectations, is a cornerstone of modern stochastic analysis. The scope of its applications, 
^ , ranging e.g. from density estimates for solutions of stochastic differential equations to 
^ ! concentration inequalities, from anticipative stochastic calculus to "Greeks" computations 
in mathematical finance, is vividly described in the three classic monographs by Malliavin 
[Tn] . Janson [9] and Nualart [20] . 

In recent years, infinite-dimensional integration by parts techniques have found another 
fertile ground for applications, that is, limit theorems and (more generally) probabilistic 
! approximations. The starting point of this active line of research is the paper [2T], where 

the authors use Malliavin calculus in order to refine some criteria for asymptotic normality 
on a fixed Wiener chaos, originally proved in [22l[24] (see also [I2j for some non-central ver- 
sion of these results). Another important step appears in [13], where integration by parts 
on Wiener space is combined with the so-called Stein's method for probabilistic approxi- 
mations (see e.g. [6l EH]), thus yielding explicit upper bounds in the normal and gamma 
approximations of the law of functionals of Gaussian fields. The techniques introduced in 
[T3] have led to several applications and generalizations, for instance: in [H] one can find 
applications to Edgeworth expansions and reversed Berry-Esseen inequalities; [18] contains 
results for multivariate normal approximations; [T] focuses on further developments in the 
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multivariate case, in relation with quasi-sure analysis; [16] deals with infinite-dimensional 
second order Poincare inequalities; in |19| . one can find new explicit expressions for the 
density of functionals of Gaussian field as well as new concentration inequalities (see also 
[3] for some applications in mathematical statistics); in [IT], the results of [I3] are com- 
bined with Lindeberg-type invariance principles in order to deduce universality results for 
homogeneous sums (these findings are further applied in [15] to random matrix theory). 

The aim of this note is to develop yet another striking application of the infinite- 
dimensional integration by parts formula of Malliavin calculus, namely the computation 
of cumulants for general functionals of a given Gaussian field. As discussed below, our 
techniques make a crucial use of a recursive formula for moments (see relation (12.21) below), 
which is the starting point of some well-known computations performed by Barbour in [2] 
in connection with the Stein's method for normal approximations. As such, the techniques 
developed in the forthcoming sections can be seen as further ramifications of the findings 
of [I3|. 

The main achievement of the present work is a recursive formula for cumulants (see 
fl4.19p ). based on a repeated use of integration by parts. Note that cumulants of order 
greater than two are not linear operators (for instance, the second cumulant coincides with 
the variance): however, our formula (I4.19P implies that cumulants of regular functionals 
of Gaussian fields can be always represented as the mathematical expectation of some 
recursively defined random variable. We shall prove in Section [5] that this implies a new 
compact representation for cumulants of random variables belonging to a fixed Wiener 
chaos. We claim that this result may replace the classic "diagram computations" adopted 
in most of the probabilistic literature (see e.g. [HEllH], as well as [23] for a general discussion 
of related combinatorial results). 

The paper is organized as follows. In Section [2] we state and prove some useful recursive 
formulae for moments. Section [3] contains basic concepts and results related to Malliavin 
calculus. Section [4] is devoted to our main statements about cumulants on Wiener space. 
Finally, in Section [5] we specialize our results to random variables contained in a fixed 
Wiener chaos. 

From now on, every random object is defined on a common suitable probability space 
(fi,^,P). 

Acknowledgements. Part of this paper has been written while we were visiting the 
Department of Statistics of Purdue University, in the occasion of the workshop "Stochastic 
Analysis at Purdue", from September 28 to October 1, 2009. We heartily thank Professor 
Frederi G. Viens for his warm hospitality and generous support. 

2 Moment expansions 

The starting point of our analysis (Proposition 12.21) is a well-known recursive relation 
involving moments and cumulants. As already discussed, this result is the seed (as well as 
a special case) of some remarkable formulae by A.D. Barbour [21 Lemma 1 and Corollary 



2 



1], providing Edgeworth-type expansions for smooth functions of random variables with 
finite moments. Since we only need Barbour's results in the special case of polynomial 
transformations, and for the sake of completeness, we choose to provide a self-contained 
presentation in this simpler setting. See also Rotar' [27] for further extensions of Barbour's 
findings. 

Definition 2.1 (Cumulants) Let X be a real- valued random variable such that < 
oo for some integer m ^ 1, and define (f)x{t) = -^(6**^), t G M, to be the characteristic 
function of X. Then, for j = 1, ...,m, the jth cumulant of X, denoted by Kj(X), is given 

by 



n,iX) = i-zy — logMt)\t=o. 



(2.1) 



For instance, ki{X) = E{X), K2{X) = E{X^) - E{Xf = Var(X), n^iX) = E{X^) 
3E{X^)E{X) + 2E{Xf, and so on. 

The following relation is exploited throughout the paper. 
Proposition 2.2 Fix m = 0, 1, 2..., and suppose that i^lXl'"^^ < oo. Then 



E{X 



m+l\ 



E 

s=0 



Ks+i{X)E{X^-^). 



Proof. By Leibniz rule, one has that 



E{X' 



-I) 



\m+l 



.s=0 



dt"^+^ 



x{t)\t=o 



d 



df^ 



-\og(j)x{t)]Mt) 



s / dt'+^ 



t=0 

(log^xW) X H) 



dt"- 



:Mt) 



(2.2) 



i=0 



with ^ equal to the identity operator. This yields the desired conclusion. 



Finally, observe that (12.21) can be rewritten as 

s=0 



E{X 



ra+l\ 



-m(m 



si 



m - s + l)E{X' 



implying (by linearity) that, for X as in Proposition l2.2l and for every polynomial / : 
of degree at most m ^ 1, 



E{xfm = E 



E 



s=0 



s\ 



d' 
dx' 



■fix) 



s=0 



df_ 

dx' 



■fix) . 
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Remark 2.3 In [21 Corollary 1], one can find sufficient conditions ensuring that the infinite 
expansion 

oo 



=0 

holds for some infinitely differentiable function / which is not necessarily a polynomial. 

3 Malliavin operators and Gaussian analysis 

We shall now present the basic elements of Gaussian analysis and Malliavin calculus that 
are used in this paper. The reader is referred to [9l [TOl [20] for any unexplained definition 
or result. 

Let ^ be a real separable Hilbert space. For any g ^ 1, let S^'^'^ be the qth tensor 
power of and denote by S^®'^ the associated gth symmetric tensor power. We write 
X = {X{h),h e 9)} to indicate an isonormal Gaussian process over S^, defined on some 
probability space P). This means that X is a centered Gaussian family, whose 

covariance is given by the relation E [X{h)X{g)] = {h, g)^. We also assume that JF = o'(X), 
that is, T is generated by X. 

For every q ^ 1, let Hg be the gth Wiener chaos of X, defined as the closed linear 
subspace of L^(f2, JF, P) generated by the family {Hq{X{h)),h e S^, \\h\\^ = 1}, where Hg 
is the gth Hermite polynomial given by 

if,(x) = (-ire4^(e-*). 

We write by convention Hq = M. For any g ^ 1, the mapping Iq^h®'^) = q\Hq{X{h)) can 
be extended to a linear isometry between the symmetric tensor product i^®^ (equipped 
with the modified norm v^lMU®?) ^ind the gth Wiener chaos Hq. For g = 0, we write 
/o(c) = c, c E M.. It is well-known that L^(f2) := L^{Q,J^, P) can be decomposed into the 
infinite orthogonal sum of the spaces Hq. It follows that any square integrable random 
variable F G L'^{Q) admits the following (Wiener-Ito) chaotic expansion 

oo 

F = Y.W,), (3.3) 

<7=0 

where /o = E[F], and the fq G S)®'^, q ^ 1, are uniquely determined by F. For every 
g ^ 0, we denote by Jq the orthogonal projection operator on the gth Wiener chaos. In 
particular, if P e L'^{Q) is as in (|3.3p . then JqF = Iq{fq) for every g ^ 0. 

Let {cfc, A; ^ 1} be a complete orthonormal system in S^. Given / G i^®^ and g G i)®'', 
for every r = 0, . . . ,pAg, the contraction of / and g of order r is the element of i^®*^^"'"'''"^''^ 
defined by 

oo 

f(g),.g= ^ {f,ei,(S)...(S)ei^)sj<8>r(g){g,ei,(S)...(g)ei^)sj<sr. (3.4) 
n,...,ir=l 
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Notice that the definition oi f ®rg does not depend on the particular choice of {cfc, k ^ 1}, 
and that f tS)r 9 is not necessarily symmetric; we denote its symmetrization by f^rQ ^ 
^0(p+cj-2r)_ Moreover, f ®og = f ® g equals the tensor product of / and g while, for p = 

f ®q9= (/,fi')s5«'J- 

It can also be shown that the following multiplication formula holds: if / G ij®^ and 
g e i:®^, then 

W)k{9) = flM^), (l)lp+,-2rmrg). (3.5) 



We now introduce some basic elements of the Malliavin calculus with respect to the 
isonormal Gaussian process X. Let S be the set of all cylindrical random variables of the 
form 

F = ^7(X(0i),...,X(0„)), (3.6) 

where n ^ 1, (7 : — > M is an infinitely differentiable function such that its partial 
derivatives have polynomial growth, and (pi E Sj, i = 1, . . . ,n. The Malliavin derivative of 
F with respect to X is the element of L'^{Q,Sj) defined as 

1=1 * 

In particular, DX{h) = h for every h E S^. By iteration, one can define the mth derivative 
D'^F, which is an element of L^(r2,i3®™), for every m ^ 2. For m ^ 1 and p ^ I, ©"^'^ 
denotes the closure of S with respect to the norm || ■ \\m,p, defined by the relation 



= E[\F\P] + Y,E{\\D^F\\l^^). 



i=l 



One also writes B>°° = flm^i flp^i D™'^- The Malliavin derivative D obeys the following 
chain rule. If (p : W ^ M is continuously differentiable with bounded partial derivatives 
and if F = (Fi, . . . , F„) is a vector of elements of then (p{F) G D^'^ and 



n 



i=l 



dxi 



Note also that a random variable F as in (13. 3p is in D^'^ if and only if 5ll'^g-^lli2(n) < ^ 
and, in this case, FdlDFUD = Xl^i ?ll-^g-^lli2(n)- If = L'^{A,£/,fi) (with /i non- 
atomic), then the derivative of a random variable F as in (13. 3p can be identified with the 
element of L^(A x Q) given by 

00 

D,F = ^gJ,„i (/,(-, x)), xgA. (3.7) 

q=l 
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We denote by 6 the adjoint of the operator D, also called the divergence operator. A 
random element u G L'^{fl,S)) belongs to the domain of 6, noted Dom5, if and only if it 
verifies \E{DF,u)sj\ ^ c„ ||-F||/^2(-q) for any F G D^'^, where Cu is a constant depending only 
on u. If M G DomS, then the random variable 6{u) is defined by the duality relationship 
(called integration by parts formula) 

E{FS{u)) = E{DF,u)s,, (3.8) 

which holds for every F G D^'^. 

The family {Pt, t ^ 0) of operators is defined through the projection operators Jg as 

oo 

Pt = Y,e-^'Jg, (3.9) 

9=0 

and is called the Ornstein-Uhlenbeck semigroup. Assume that the process X', which stands 
for an independent copy of X, is such that X and X' are defined on the product probability 
space {Q x ^ (g) P x P'). Given a random variable F G D^'^, we can regard it as 
a measurable mapping from to M, determined P o X~^-almost surely. Then, for any 
t > 0, we have the so-called Mehler's formula: 



PtF = E'{F{e-'X + Vl - e-2*X')), (3.10) 

where E' denotes the mathematical expectation with respect to the probability P'. By 
means of this formula, it is immediate to prove that Pt is a contraction operator on L^(f2), 
for all p ^ 1 . 

The operator L is defined as L = X]^o ~1'^qi it can be proven to be the infinitesimal 
generator of the Ornstein-Uhlenbeck semigroup {Pt)t^o- The domain of L is 



,2,2 



DomL = {F e L^{n) : UgHl^n) < oo} = ©' 

q=l 

There is an important relation between the operators D, S and L. A random variable F 
belongs to D^-^ if and only if F G Dom {6D) (i.e. F G D^'^ and DF G Dom5) and, in this 

6DF = -LF. (3.11) 

For any F G L'^{^), we define L~^F = The operator is called the 

pseudo-inverse of L. Indeed, for any F G L'^{^), we have that L~^F G DomL = D^'^, and 

LL-^F = F - E{F). (3.12) 

We now present two useful lemmas, that we will need throughout the sequel. The 
first statement exploits the two fundamental relations (13. lip and f l3.12p . Note that these 
relations have been extensively applied in [I3|, in the context of the normal approximation 
of functionals of Gaussian fields by means of Stein's method. 
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Lemma 3.1 Suppose that F G D^'^ and G G L'^iyt). Then, L~^G G ©^'^ and we have: 

E{FG) = E{F)E{G) + E{{DF, -DL-^G)^). (3.13) 

Proof. By (fXTTH and (jmi) . 

E{FG) - E{F)E{G) = E{F{G - E{G))) = E{F x LL'^G) = E{F6{-DL-^G)), 

and the result is obtained by using the integration by parts formula (|3.8I) . ■ 

Remark 3.2 Observe that {DF, —DL'^G) is not necessarily square-integrable, albeit it 
is by construction in L^{Q). On the other hand, we have 

POO 

{DF,-DL^^G)^= / e~'{DF,PtDG)^dt, (3.14) 
Jo 

see indeed identity (3.46) in [T3] . 

The next result is a consequence of the previous Lemma ISTTl 

Proposition 3.3 Ftx p^2, and assume that F G L^P-^{n) n P^'^. T/ten, F^ G D^'^ and 
= pFP-^DF. Moreover, for any G G L'^i^), 

E{FPG) = E{FP)E{G) +pE{FP-\DF,-DL-^G)s^). (3.15) 

4 A recursive formula for cumulants 

The aim of this section is to deduce from formula (12.21) a recursive relation for cumulants 
of sufficiently regular functionals of the isonormal process X. To do this, we need to 
(recursively) introduce some further notation. 

Definition 4.1 Let F G D^'^. We define ro(F) = F and ri(F) = {DF, -DL^^F) If, 
for j ^ 1, the random variable Tj{F) is a well-defined element of L'^{fi), we set Tj+i{F) = 
{DF, —DL^^Vj{F))f^. Observe that the definition of rj+i(F) is well given, since (as already 
observed in general) the square-integrability of Tj{F) implies that L~^Vj{F) G DomL = 

The following statement provides sufficient conditions on F, ensuring that the random 
variable Tj{F) is well defined. 

Lemma 4.2 1. Fix an integers j ^ 1, and let F,G e D^'^' . Then {DF, ~DL~^G)sj G 
D-'~^'2^ , where we set by convention (but consistently!) ID)^'^ = L^{fl). 

2. Fix an integer j ^ 1, and let F G W'"^^ . Then, for all k = 1, . . . ,j, we have that Tk{F) 
is a well-defined element ofW"^'"^^ ; in particular, one has that Tk{F) G L^{VL) and 
that the quantity £'[rfc(F)] is well-defined and finite. 
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3. If F G ( in particular if F equals a finite sum of multiple integrals ), then Tj (F) G 
poo every j ^ 0. 

Proof. Without loss of generality, we assume during all the proof that 9) has the form 
L^(A, /i), where {A,^/) is a measurable space and // is a a-finite measure with no 
atoms. 

1. Let k E {0, . . . , j — 1}. Using Leibniz rule for D (see e.g. [20l Exercice L2.13]), one 
has that 



- D^{DF,-DL~^G)^ = I DaF DaL-^G fi{da) 

J A 



(4.16) 



J2 (^) J^D^~\DaF)^D\DaL-'G)n{da), 



with ® the usual symmetric tensor product. Note that, to deduce ( ]4.16p . it is sufficient to 
consider random variables F, G that have the form (13. 6p (for which the formula is evident, 
since it basically boils down to the original Leibniz rule for differential calculus), and then 
to apply a standard approximation argument. From (I4.16p . one therefore deduces that 

\\D''{DF,-DL-^G)f,hm 



1=0 
k 



D^^\DaF)®D\DaL~^G) ^i{da) 



1=0 
k 

^ E 

1=0 

k 



^ |JD''-\DaF)\\^^,.^4D\DaL-'G)\\^^l^I{do 



E 

1=0 



' ljD>^-^iD,F)\\l^,,_,f^{da)J 1^^^^^ 



^«i{k-l+l) 



(4.17) 



By mimicking the arguments used in the proof of Proposition 3.1 in [16] (see also (13.141) ). 
it is possible to prove that 



-D'+^L-'G 
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Consequently, for any real p ^ 1, 



{I 






1 












' F 



i)]dt 



(4.18) 



where, to get the last inequality, we used the contraction property of Pt on L'^^fl). Fi- 
nally, by combining (14.171 ) with (14.181) through the Cauchy-Schwarz inequality on the one 
hand, and by the assumptions on E and G on the other hand, we immediately get that 
\\D''{DE,-DL-'^G)^\\^m belongs to L'^'~\n) for all k = 0, - 1, yielding the an- 
nounced result. 

2. Fix j ^ 1 and E G The proof is achieved by recursion on k. For k = 1, the 

desired property is true, due to Point 1 applied to G = F. Now, assume that the desired 
property is true for k j — 1), and let us prove that it also holds for k + 1. Indeed, we 



have that Tk+iiE) = {DE, -DL-^Vk{E)) f,, that E e W'"^' C 



(assumption on E) 



and that Tk{E) G 



(recurrence assumption). Point 1 yields the desired conclusion. 



3. The proof is immediately obtained by a repeated application of Point 2. 



□ 



We will now provide a new characterization of cumulants for functionals of Gaussian 
processes: it is the fundamental tool yielding the main results of the paper. Note that, due 
to Lemma [42] (Point 2), the following statement applies in particular to random variables 



m 



(m ^ 2). 



Theorem 4.3 Fix an integer m ^ 2, and suppose that: (i) the random variable E is an 
element of L'^^~'^{yt) fl D^'^, (ii) for every j ^ m — 1, the random variable Tj{E) is in 
L'^{Q). Then, for every s ^ m, 



Ks+i{E) = s\E[T,iE)]. 



(4.19) 



Proof. The proof is achieved by recursion on s. First observe that ni{E) = E{E) = 
E[rQ{E)], so that the claim is proved for s = 0. Now suppose that (14.190 holds for every 
s = 0, /, where I ^ m — 1. According to (12.21) . we have that 



E{E^^') = Y,Qns^,{E)E{E^~^) 



f^i+i{E). 



(4.20) 



On the other hand, by applying (13.15P to the case p = I and G = F, we deduce that 

E{E^+^) = E{E^)E{E) + IE{E^~^DE, -DL~^E)^) = E{E^)ki{E) + lE{E^~^Ti{E)). 
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By the recurrence assumption, and by applying again (I3.15P to the case p = I — 1 and 
G = ri(F), one deduces therefore that 

= E{F^)ki{F) + lE{F^~^Vi{F)) 
= t^r{F)E{F') + K2{F)IE{F'-^) + /(/ - l)E{F'-^T2{F)). 

Iterating this procedure yields eventually 



E{F'+') = J2 (^)Ks+iiF)EiF^-^) + l\E[TiiF)], 



SO that one deduces from fl4.20p that relation (I4.19p holds for s = / + 1. The proof is 
concluded. ■ 

Remark 4.4 Suppose that F = Jg(/), where g ^ 2 and / G i^®''. Then, L'^F = -q~'^F, 
and consequently 

E[T,{F)] = E[{DF,-DL~'Ts^i{F))^]=E[FT,_i{F)] (4.21) 
= E[{-DL~'F,Dr,.,{F))^] = ^E[{DF,Drs^i{F))s,]. 

Several applications of formula (14.191) (and, implicitly, of f l4.2ip ) are detailed in the 
next section. 



5 Cumulants on Wiener chaos 
5.1 General statement 

We now focus on the computation of cumulants associated to random variables belonging 
to a fixed Wiener chaos, that is, random variables having the form of a multiple Wiener-Ito 
integral. Our main findings are collected in the following statement, providing a quite com- 
pact representation for cumulants associated with multiple integrals of arbitrary orders. In 
the forthcoming Section 15.21 we will compare our results with the classic diagram formulae 
that are customarily used in the probabilistic literature. 

Theorem 5.1 Let q ^ 2, and assume that F = Iq{f), where f G S)'^'^. Denote by Hs{F), 
s ^ 1, the cumulants of F. We have i^i{F) = 0, K2{F) = glU/Hl®, and, for every s ^ 3, 

K,{F) = q\{s-l)\ ^ c,(ri, . . . , r,_2)((...((/®. J)®. J) . . . ®.,_3/)®r,. J, />^«„ (5.22) 
where the sum ^ runs over all collections of integers ri, . . . ,rs_2 such that: 
(i) 1 ^ ri, . . . ,r,„2 ^ q; 
(n) ri + ...+r,_2 = ^; 
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(iii) ri < g, ri + r2 < f , ri + ... + r^_3 < 

(iv) r2 ^ 2g - 2ri, . . ., r^.g ^ (s - 2)g - 2ri - . . . - 2r,_3; 

and where the combinatorial constants Cg(ri, . . . ,rs_2) are recursively defined by the rela- 



tions 

in I \ ^ 



q-l 
r — 1 



and, for a ^ 2, 



I N / , ag - 2ri - . . . - 2ra_i - 1\ /g - 1 , . 

Cq(ri, . . . ,r„) = g(r„ - 1)! , Cg(ri,...,r„_i, 



Remark 5.2 1. If sg is odd, then ^^(F) = 0, see indeed condition {ii). 

2. If g = 2 and F = hif), f ^ then the only possible integers ri, . . . , rs_2 verifying 
(i) — (if) in the previous statement are ri = . . . = rs_2 = 1. On the other hand, 
we immediately compute that C2(l) = 2, 02(1, 1) = 4, 02(1, 1,1) = 8, and so on. 
Therefore, 

Ks{F) = r~\s - 1)!((...(/ ®i /) . . . /) ®i /, />^«„ 

and we recover the classical expression of the cumulants of a double integral (as used 
e.g. in [7| or [HI). 

3. If g ^ 2 and F = /g(/), / G i^®^ then (15:221) for s = 4 reads 

9-1 

«:4(/,(/)) = 6g!5^c,(r,g-r)((/®J)®,_J,/)^^, 

r=l 

= ^ E (^) '(2? - 20!((/® J) ®,~r /, />^«, 

r=l ^ ^ 

r=l ^ ^ 

= -E^^!'(^)'(2g-2r)!||/^,/|||,(,,_,.), (5.23) 

and we recover the expression for ^^{F) deduced in [T71 Section 3.1] by a different 
route. Formula f l5.23p should be compared with the following identity, first estab- 
lished in [22l p. 183]: for every g ^ 2 and every / G 



/2g-2r\ ~ ,,2 
g — r y " 
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Note that it is in principle much more difficult to deal with (15.241) . since it involves 
both symmetrized and non-symmetrized contractions. 

Proof of Theorem \5.1\ Let us first show that the following formula (15.251) is in order: for 
any s ^ 2, 

q [(s-l)g-2ri~...-2rs-2]A(j 
ri=l rs_i=l 

X I sq—2ri — ...—2rs-i_ ((•••(/®n/)®rj)... J). 

(5.25) 

We shall prove (15.251) by induction, assuming without loss of generality that has the 
form L^(A, /i), where (A, s^) is a measurable space and /i is a cr-finite measure without 
atoms. When s = 2, identity (I5.25P simply reads 

9 

ri(F) = Y,^q{r)hq-2rU®rf). (5.26) 
r=l 

Let us prove (I5.26p by means of the multiplication formula (13.51) , see also [2T] for similar 
computations. We have 

ri(F) = {DF,-DL-^F)^=-\\DF\\l = q f I,^^{f{-,a)fii{da) 

Q J A 

= hq~2~2r(^j /(", a)®r/(-, a)/i(rfa)^ 

= q^rli J /2<?-2-2r (/®r+l/) 

= gE(^-i)!(^:!)V2r(/®./)' 

r=l ^ ^ 



thus yielding (15.261) . Assume now that (I5.25p holds for Ts_i{F), and let us prove that 
it continues to hold for Ts{F). We have, still using the multiplication formula ( 13.51) and 
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proceeding as above, 

q [{s~l)q-2ri-...-2ra-2\M 

r-i=l r-s_i=l 

Xl{ri + ...+r,_i<f }(Vl(/)>^^g-2n-...-2rwi-l((...(/®n/)®r2/) .../)®r,_i/)>^ 

q [(s-l)(y-2ri-...-2rs_2]Ag [sq~2ri~...~2rs-i]Aq 

= Y Cg(ri,...,r,„i) X g(r, - 1)! 

ri=l rs-.i=l rs=l 

/sg-2ri-...-2r,_i-l\ 

Xl{r-i + ...+r-,_i<f }A^+l)g-2ri-...-2r, ((• • • (/®n /)®r2 /) • • • /)®rj), 

which is the desired formula for Ts{F). The proof of (I5.25p for all s ^ 1 is thus finished. 
Now, let us take the expectation on both sides of (j5.25l) . We get 

Ks{F) = (s-l)!E[r,_i(F)] 

q [(s-l)Q-2ri-...-2rs_2]Ag 

= Cg(ri,...,r,_i)l|,, <,}... 

ri=l rs_i=l 

Observe that, if 2ri + . . . + 2rs_i = sg and Vs-i ^ (s — l)g — 2ri — ... — 2rs_2 then 
2rcj_i = g + (s — l)q — 2ri — ... — 2rs__2 ^ Q + rs~i, so that r^^i ^ g. Therefore, 

q [{s-2)g-2ri-...-2r^_3]Ag 



ri=l rs-2=l 

xl 



which is the announced result, since Cg{ri, . . . , rs_2, q) = q\Cg{ri, . . . , rs-2)- □ 



5.2 Combinatorial expression of cumulants 

We now provide a classic combinatorial representation of cumulants of the type Ks{F), in 
the case where: (i) s ^ 2, (ii) q ^ 2, (iii) sq is even, (iv) F = Iq{f) (with / e ^®^) and 
(v) S) = L'^{A,s^,fi) is a non-atomic measure space. Assumptions (i)-(v) will be in order 
throughout this section. As explained e.g. in |23| . one can equivalently express cumulants 
of chaotic random variables by using diagrams or graphs: here, we choose to adopt the 
(somewhat simpler) representation in terms of graphs. See [23l Section 4] for an explicit 
connection between graphs and cumulants; see [H] for some striking application of graph 
counting to the computation of cumulants of non-linear functionals of spherical Gaussian 
fields; see [H [5], [8] for classic examples of how to use diagram enumeration to deduce CLTs 
for Gaussian subordinated fields. 
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Definition 5.3 For s ^ 2, consider the set [s] = {1, s} of the first s integers. For g ^ 2, 
we denote by K{s, q) the class of all non-oriented graphs 7 over [s] satisfying the following 
properties: 

- 7 does not contain edges of the type (j, j), j = 1, s, that is, no edges of 7 connect 
a vertex with itself. One can equivalently say that 7 "has no loops". 

- Multiple edges are allowed, that is, an edge can appear k ^ 2 times into 7; in this 
case, we say that k is the multiplicity of the edge. Also, if i,j are connected by an 
edge of multiplicity k, we say that i and j are connected k times. 

- Every vertex appears in exactly q edges (counting multiplicities). 

- 7 is connected. 

If sq is odd, then K{s, q) is empty. If sq is even, then each 7 G K{s, q) contains ex- 
actly sq/2 edges (counting multiplicities). For instance: an element of -ft"(3,2) is 7 = 
{(1, 2), (2, 3), (3, 1)}; an element of K(3, 4) is 7 = {(1, 2), (1, 2), (2, 3), (2, 3), (1, 3), (1, 3)} 
(note that each edge has multiplicity 2). 

Definition 5.4 Given g, s ^ 2 such that sq is even, and 7 G K{s, g), we define the constant 
u'(7) as follows. 

(a) For every j = 1, s, consider a generic vector L{j) = (/(j, 1), q)) of g distinct 
objects. Write {L{j)} for the set of the components of L{j). 

(b) Starting from 7, build a matching 772(7) over L := IJj=i s{-^(i)} follows. Enu- 
merate the vertices Vi, ..■,Vsq/2 of 7. If links ii and ji and has multiplicity ki, then 
pick ki elements of L{ii) and ki elements of L{ii) and build a matching between the 
two /ci-sets. If f2 links 12 and j2 and has multiplicity ^2, then pick k2 elements of 
L{i2) (among those not already chosen at the previous step, whenever 12 equals ii or 
ji) and k2 elements of L{i2) (among those not already chosen at the previous step, 
whenever ^2 equals ii or ji) and build a matching between the two /c2-sets. Repeat 
the operation up to the step sq/2. 

(c) Define &q as the group of all permutations of [q]. For every a G ©g, define the vector 
LaU) = 1), ...Ja{j,q)), where la{j,p) = /(jXp)), P = 1, 

(d) Define as the sth product group of &q, that is, &q^^ is the collection of all s- 
vectors of the type a = (cti, 0"^), where aj G &q, j = l,...,s, endowed with the 
usual product group structure. 

(e) For every a = (ai, ...,0"s) G build a new matching rricr^'y) over L by repeating 
the same operation as at point (b), with the vector Lfj^{j) replacing L{j) for every 
j = l,...,s. 



14 



(f) Define an equivalence relation over &[^^ by writing a ~^ vr whenever mo-(7) = 771,^(7). 

(s) (s) 

Let &q / be the quotient of &q with respect to 

(g) Define ^(7) to be the cardinality of ©g^'V 

For instance, one can prove that 1^(7) = 2"*"^ for every s ^ 2 and every 7 G 2). Also, 
^(7) = q\ for every q ^ 2 and every 7 G -ft'(2, g). 

Definition 5.5 For g ^ 2, let / G i^®'', that is, / is a symmetric element of L^(A'', ^Z'^, /i'^). 
Fix 7 G K{s, q), where s ^ 2 and sq is even. Starting from / and 7, we define a function 

(ai, asg/2) ^ f-yiai, 
of sq/2 variables, as follows: 

(i) juxtapose s copies of /, and 

(ii) if the vertices j and / are linked by r edges, then identify r variables in the argument 
of the jth copy of / with r variables in the argument of the Ith copy. By symmetry, 
the position of the identified r variables is immaterial. Also, by connectedness, one 
has necessarily r < q. 

The resulting function is a (not necessarily symmetric) element of 

For instance, if 7 = {(1,2), (2,3), (3,1)} G i^(3,2), then f^{x,y,z) = f{x,y)f{y,z)f{z,x). 
If 7 = {(1, 2), (1, 2), (2, 3), (2, 3), (1, 3), (1, 3)} G K{3, 4), then 

f^{t, u, V, X, y, z) = fit, u, V, x)f{t, u, y, z)f{y, z, v, x). 

We now turn to the main statement of this section, relating graphs and cumulants. 
The first part is classic (see e.g. p3| for a proof), and shows how to use the functions 

to compute the cumulants of the random variable F = Iq{f)- The second part of the 
statement makes use of (15.221) . and shows indeed that Theorem 15.11 implicitly provides a 
compact representation of well-known combinatorial expressions. 

Proposition 5.6 Let q ^ 2 and assume that F = Iq{f), where f G i^®''. Assume the 
integer s ^ 2 is such that sq is even. Then, 

Ks{F) = w{-f) f^{ai, ...,asq/2)f^{dai) ■ ■ ■ i2{dasq/2)- (5.27) 

l<^K{s,q) 

As a consequence, by using ( f5. 2S\) . one deduces the identity 

(s - ly.ql J2 Cqin, . . . , r,_2)((...(/®.,/)®.,/) . . . /)®.._ J, />^«, (5.28) 

= w(7) / f^{ai, ...,asq/2)lJ'{dai) ■ ■ ■ j2{dasq/2)- (5.29) 

7GX(s,g) 

where J2 funs over all collections of integers ri, . . . , rs-2 verifying the conditions pinpointed 
in the statement of Theorem \5.1[ 
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Remark 5.7 Being based on a sum over the whole set K{s, g), formula (I5.27P is of course 
more compact than fl5.22l) . However, since it does not contain any hint about how one 
should enumerate the elements of K{s,q), expression (15.271) is much harder to compute 
and asses. On the other hand, (15.221) is based on recursive relations and inner products 
of contractions, so that one could in principle compute Ks(-F) by implementing an explicit 
algorithm. 

5.3 CLTs on Wiener chaos 

We conclude the paper by providing a new proof (based on our new formula (j5.22l) ) of the 
following result, first proved in [22] and yielding a necessary and sufficient condition for 
CLTs on a fixed chaos. 

Theorem 5.8 (See |22| ) Fix an integer q ^ 2, and let (-F„)n^i be a sequence of the form 
Fn = Igifn), with fn E ^i®"^ such that E[F^] = q\\\fn\\lm =J for all n ^ 1. Then, 
as n ^ oo, we have F^ N[0, 1) in law if and only if \\fn^rfn\\s}«>(2q-2r) for all 
r = 1, . . . ,g - 1. 

Proof. Observe that = and K2(i^„) = 1. For s ^ 3, we consider the 

expression (j5.22l) . Let ri,...,rs_2 be some integers such that {i)-{iv) in Theorem 15.11 
are satisfied. Using Cauchy-Schwarz inequality and then successively \\g<S)rh\\^<g>{p+q-2r) ^ 
\\g ®r h\\ ^ts{p+q~2r) ^ ||fi'|U8'p||^|U«9 whenever g G i^®^, h G and r = 1, . . . ,p A g, we get 
that 

I {{■■■{fn®rifn)®r2fn) ■ ■ ■ fn)®rs-2fn^ fn)sj'^i \ 

^ ||(---(/n®ri/n)®r2/ri) • • • /ra) ®r,_2 /n || /n || S5®9 

ii/niir.', 

= (?!)^"^ ||/n®ri/n|U®(29-2ri). (5.30) 

Assume now that || /n®r /nlU®(29-2r) — for all r = 1, . . . , g — 1, and fix an integer s ^ 3. 
By combining (I5.22p with (15.301) . we get that Ks(-F„) ^ as n — > oo. Hence, applying the 
method of moments or cumulants, we get that F„ — > A^(0, 1) in law. Conversely, assume 
that Fn —>■ N N{0, 1) in law. Since the sequence (F„) lives inside the gth chaos, and 
because E[F^] = 1 for all n, we have that, for every p ^ 1, sup„>j^ -EHF^p] < cxo (see e.g. 
Janson [9l Ch. V]). This implies immediately that n^iFr,) = E[F^] - 3 ^ E[N^] -3 = 0. 
Hence, identity (|5.23l) allows to conclude that \\fn®rfn\\f)<s{2q-2r) — * for all r = 1, . . . , g — 1. 

□ 
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